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SUMMARY 

This report presents the latest modifications made to the computer code TURBAN, which 
does a preliminary sizing analysis for axial-flow turbines. The TURBAN analysis is based on 
mean-diameter flow characteristics. Program input includes flow, speed, and power or 
pressure ratio. The output presents annulus dimensions, diagram velocities and angles, and 
efficiencies. Options are provided for varying stage number, mean diameter, reaction, 
loading, diagram type, and/or work split. 

Modifications were recently made to TURBAN to satisfy user needs and convenience. 
Turbine cooling-air flows and temperature now can be accounted for along with an 
associated efficiency decrement. Alternative input options have been added for defining the 
velocity diagram by stage reaction, for setting the mean diameter by stage loading, and for 
arbitrarily specifying stage work split. The Reynolds number dependency for the loss model 
was weakened, and an internal calculation of air viscosity was added as a default. The 
analytical modeling for these modifications are presented herein. 

This report also serves as an updated users manual for the modified TURBAN code. 
Program input and output are described, and sample cases are included for illustration. 

INTRODUCTION 

Preliminary studies of gas turbine systems require many repetitive calculations of 
geometry and design-point performance for all the components. For this type of screening 
analysis, rapid approximate procedures, rather than time-consuming rigorous procedures, 
are sufficient to yield the desired component overall geometry and performance 
characteristics. One such analysis code, named TURBAN, for axial-flow turbines was first 
presented in reference 1 more than 20 years ago. An updated version of this code with 
numerous modeling improvements was presented in reference 2, and blade geometry 
modeling subsequently added to the code was reported in reference 3. 

Recent use of TURBAN for aircraft propulsion system studies resulted in the desire for 
additional capabilities that were then added to the code. Turbine cooling now can be includ- 
ed in the analysis. New alternative input options allow direct specification of stage reaction, 
stage loading, and stage work distribution. The Reynolds number dependency was modified 
to provide improved loss modeling, and an internal calculation of air viscosity was added for 
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convenience. These modifications require the use of additional input information. 

This report presents the analytical modeling associated with the new capabilities added to 
the TURBAN code. It also serves as an updated users manual for the code. Program input 
and output are described. Sample cases are included for illustration. 

SYMBOLS 


c p heat capacity, J/kg-K; Btu/lb-°R 

D mean diameter, m; ft 

g dimensional constant, 1; 32.17 lbm-ft/lbf-sec 2 

Ah specific work, J/kg; Btu/ib 

J dimensional constant, 1 ; ft-lb/Btu 

j stage number for change in meanline slope 

m number of stages for meanline diameter variation 

N rotative speed, rad/sec; rpm 

n number of turbine stages 

P shaft power, W; hp 

p total pressure, N/m 2 ; lb/ft 2 

R stage reaction 

Re Reynolds number 

T total temperature, K; °R 

U blade speed, m/sec; ft/sec 

V gas absolute velocity, m/sec; ft/sec 

W gas relative velocity, m/sec; ft/sec 

w mass flow rate, kg/sec; Ib/sec 

x ratio of stage work to turbine work 

y ratio of blade-row coolant flow to blade-row inlet flow 

5 stage-efficiency specific loss 

Sr relative roughness 

il total efficiency, overall or stage 

p. viscosity, N-sec/m 2 ; Ib/sec-ft 

\y stage work factor 

Subscripts: 
c coolant 

ex turbine exit 

i index for stage number 

in turbine inlet 

n last stage 

p primary (turbine inlet) 

st stator 
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ro rotor 

rp rotor primary 

u tangential component 

1 first stage or stator exit 

2 rotor exit 

Superscript: 

y specific heat ratio 

* corrected for coolant 


ANALYTICAL MODELING 

The analytical modeling for the enhanced capabilities are presented in this section. The 
added models used for the turbine cooling, velocity diagram, and loss model calculations are 
discussed. 


Turbine Cooling 

The model used to account for turbine cooling is based on mixing the primary flow and all 
the cooling flow (assuming constant heat capacity) at the turbine inlet to obtain a corrected 
turbine inlet temperature. 


T in = (WpTin + Sw c>i T c ) / (w p + Ew c .|) ( 1 ) 

Both the array of blade-row cooling flows, w ci , and the coolant temperature, T c , are required 
as code input. This methodology is consistent with the turbine thermodynamic efficiency 
definition based on cooling air pressure being equal to turbine inlet pressure. 

Tlth = P / {(WpT in + Iw c>i Tc)Cp[1 - (Pex/Pin) (r1)/Y ]} ( 2 ) 

The turbine efficiency provided by the loss model is for an uncooled turbine. To determine 
cooled turbine thermodynamic efficiency, reference 4 uses assigned values of 
stage-efficiency specific loss, which is defined as 


6 = Ati / (r| unc y) ( 3 ) 

for stator and rotor to determine the stage-efficiency loss due to cooling. The cooled stage 
efficiency is then obtained from the uncooled stage efficiency as 


titf/Hunc — 1 ‘ y S tS s t " yro^ro ( 4 ) 

Values for stage-efficiency specific loss, which are code inputs for stator and rotor, are given 
in reference 4 for various cooling configurations. These values are based on limited data and 
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20-year-old technology. Values relevant to current technology do not appear to be generally 
available. 

Cooled turbine efficiency is often alternately expressed as rotor primary efficiency, which 
is defined as 


% = p / {(WpT ln + W c1 T c )Cp[1 - (Pe/p»,) lr, ^J} ( 5 ) 

Rotor primary efficiency can be obtained from thermodynamic efficiency by combining 
equations (2) and (5) for the same output power. 

hrp = •nth(w p T in + Xw Cii T c ) / (WpT in + w c1 T c ) ( 6 ) 

Velocity Diagram Options 

Options have been added to the code that influence the velocity diagrams. These include 
specifying the stage reaction to define the velocity diagram, specifying the stage work factor 
to define the mean diameter, and arbitrarily specifying the stage work split. 

Stag e reaction input.- The definition of stage reaction, assuming constant blade speed 
across a rotor, is 


R = (W 2 2 -W 1 2 ) / ( V ! 2 -V 2 2 +W 2 2 -W 1 2 ) 


(7) 


and the definition of stage work factor is 


\j/ = AV u /U (8) 

With the assumption of constant axial velocity across the stage, equations (7) and (8) can be 
combined with the velocity diagram equations to relate stage swirl split to stage reaction and 
work factor. 


V u ,i/AV u = (1 -R)/\|/ + 0.5 (9) 

Therefore, an input value of stage reaction, R, serves to define the velocity diagram since the 
stage work factor is also known from the input. 

Stage work-factor input-- A stage-average work factor can be defined based on the mean 
squared blade speed. 


\j/ = gJIAhj / EU 2 (10) 

By inputting this stage work factor, the stage mean blade-speed summation is computed 
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from equation (10) and the stage mean diameter summation from 

ED) 2 = (720 / jtN) 2 EU, 2 (11) 

For a single-stage turbine or a constant mean-diameter multistage turbine, the single 
value of mean diameter, Dj, can be calculated directly from equation (11). With a varying 
mean diameter, the mean diameter summation is expressed 

ED 2 = D 1 2 E(D i /D 1 ) 2 (12) 

Since the mean diameter variation is linear with stage number (ref. 2), the summation can be 
expressed using arithmetic progression summation formulas in terms of the exit to inlet 
diameter ratio, Dp/Di, which must be input, and the number of stages, m, over which the 
diameter variation occurs. 

E(D, / DO 2 = /(Dn/D^m) = miDp/Di + (Dp/D, - 1) 2 (2m - 1) / [6(m - 1)]} ( 13 ) 

Corresponding to the three options (ref. 2) available for stage mean-diameter variation, the 
diameter-ratio summation for equation (12), using the function notation defined by 
equation (13) becomes: 

(1) Linear variation between first and last stages 

E(Dj / Di) 2 =/(Dn/D 1p n) (14) 

(2) Constant from first stage to j ,h stage and then linear to last stage 

S(D, / D,) 2 = [j - 1 + /( Dn/D^n+l-j)] ( 15 ) 

(3) Linear from first stage to jth stage and then constant to last stage 

E(D, / D-i) 2 = [(n-j)(D n / D^ 2 + /(Dn/D^j)] (16) 

The inlet diameter, D 1f can now be obtained from equation (12) using ED 2 computed from 
equation (1 1) and E(D| / D^ 2 computed from equation (14), (15), or (16) as appropriate. The 
exit diameter, D n , is then found from the input diameter ratio Dp/Dv 

Stage work-sp lit input.- Previously, the basic assumption of constant stage work factor 
resulted in stage work split being determined uniquely by the stage diameter variation. 

Ah, / EAhj = U, 2 / EU, 2 = D, 2 / ED, 2 (17) 


An option is now available for the direct specification of an arbitrary work split 
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Ahj/IAhj = Xj 


(18) 


where the X| are input. As a result, stage work factor is not constant for this option 

\j/i = gJxiEAhi / U, 2 (19) 


and equation (18) replaces equation (17) as required in the analysis of reference 2. 

Loss Model 

The dependency of loss on Reynolds number was based (ref. 2) on flow in a smooth tube. 

Loss « Re 2 (20) 

Within the limited Reynolds number variation of the turbine database used for loss model 
calibration, this model appeared satisfactory. Subsequent studies of turbines with larger 
variations in Reynolds number indicated that this dependency was too strong. Therefore, it 
was replaced by an implicit approximation (ref. 5) of the Karman-Prandtl equation 

Loss - {log[6.9 / Re + (e r / 7.4) 111 ]}' 2 ( 21 ) 

A relative roughness, e r , of 0.0002 was used for equation (19). 

For user convenience, an internal calculation of viscosity for air using the Sutherland 
equation was added as a default option. 

H = 7.238 x 10' 7 T 1 5 / (T + 199) ( 22 ) 

DESCRIPTION OF INPUT AND OUTPUT 

This section serves as an updated users manual by presenting a detailed description of 
the program input and normal output. The error messages are as described in reference 2. 
Included in the input and output sections are example cases illustrating the use of the 
program with the new options. 


Input 

The program input, a sample of which is presented in table I, consists of a title record and 
the required physical data and option indicators in NAMELIST form. The title, which is print- 
ed as a heading on the output listing, can contain up to 77 characters located anywhere in 
columns 2 through 78 on the title record. A title, even if it is left blank, must be the first record 
of the input data. Additional titles can be used to identify different cases in the same data file. 
This is done by placing a title in front of the data for the particular case and using the option 
indicator ITIT as subsequently described. 
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The physical data and option indicators are input in data records having the NAMELIST 
name INPUT. The variables and indicators that comprise INPUT and the proper units are as 
follows. These must be input for all cases except where otherwise indicated. Either SI units 
or U.S. customary units may be used. 

PTIN inlet total pressure, N/cm 2 ; lb/in. 2 

TTIN inlet total temperature, K; °R 

MU gas viscosity, N-sec/m 2 ; Ib/sec-ft. 

<0.0 - internal computation of viscosity for air 
>0.0 - value of viscosity 

R gas constant, J/kg-K; ft-lbf/lbm-°R 

GAM specific heat ratio 

DIN inlet diameter - hub or mean or tip value as specified by IDIAM=1 -3, cm; in. 

- relative (to exit) value if IDIAM=4 

DEX exit diameter - hub or mean or tip value as specified by IDIAM=1-3, cm; in. 

- relative (to inlet) value if IDIAM=4 

RREX exit radius ratio; RREX may be omitted when IDIAM=2 or 4 and IALPH=0; RREX 
is used as a first trial when IALPH=0 and IDIAM=1 or 3 

RPM rotative speed, rad/sec; rpm 

POW shaft power - omit when IPR=1 , kW; hp 

W mass flow rate, kg/sec;lb/sec 

ALPHA stator exit angle from axial direction; ALPHA is used as first trial value when 
IALPH=1 , deg 

ALPHAO turbine inlet flow angle from axial direction; input only when KALPH0=2, deg 
VU1DVU ratio of rotor inlet swirl to total change in swirl; input only when IVD=5 
REACT stage reaction; input only when IVD=6 
WF stage work factor; input only when IDIAM=4 
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Xl(l) ratio of stage work to total work; Xl(1 )>0.0 triggers this option, which requires 

1=1,NMIN that NMIN=NMAX 

KLOSS turbine loss coefficient; a value of 0.3 is recommended in the absence of 
additional information 

NMIN minimum number of stages for which the calculations are performed 

NMAX maximum number of stages for which the calculations are performed; results are 
obtained for all stage numbers between NMIN and NMAX 

NMID stage number at which the meanline changes slope; input only when IMID=1 

E squared ratio of stage-exit to stage-average meridional velocities 

PRTS turbine inlet-total to exit-static pressure ratio; input only when IPR=1 

WCOWP(I) ratio of blade-row coolant flow to turbine-inlet flow; (default=NMIN*0.0) 

1=1 .NMIN 

TCOTP ratio of coolant temperature to turbine inlet temperature; 

DELS stage-efficiency specific loss for stator cooling; (default=0.1 5) 

DELR stage-efficiency specfic loss for rotor cooling; (default=0.30) 

IALPH indicates whether stator exit angle or turbine exit radius ratio is specified: 

= 0 - turbine is designed for specified ALPHA 
= 1 - turbine is designed for specified RREX 

IDIAM indicates whether input diameters are absolute hub, mean, or tip values or a 
relative mean value: 

= 1 - input diameters are hub values 
= 2 - input diameters are mean values 
= 3 - input diameters are tip values 
= 4 - input diameters are relative mean values 

IVD indicates type of velocity diagram used: 

= 1 - symmetrical diagrams 
= 2 - zero exit swirl diagrams 
= 3 - impulse diagrams 
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= 4 - zero exit swirl diagrams if \j/ < 2.0 and impulse diagrams if \|/ > 2.0 
= 5 - ratio of rotor-inlet swirl to total change in swirl is input as VU1 DVU 
= 6 - stage reaction is input as REACT 

ITIT indicates use of titles in addition to that required as first line of data file: 

= 0 - no title precedes next case 

= 1 - title line precedes next case; must be input for each additional title 
because ITIT is set to zero after each title is read 

IEV indicates use of exit vanes: 

= 0 - no exit vanes 

= 1 - exit vanes are used to turn turbine exit flow to axial direction 

IPR indicates whether shaft power or pressure ratio is specified: 

= 0 - shaft power is input 

= 1 - turbine inlet-total to exit-static pressure ratio is input 

IU indicates type of units used for input and output: 

= 1 - SI units 
= 2 - U.S.customary units 

KALPHO indicates turbine-inlet flow angle option: 

= 0 - turbine-inlet flow is axial (default) 

= 1 - turbine-inlet flow angle equals stage-exit flow angle 
= 2 - turbine-inlet flow angle is input as ALPHAO 

IAR indicates blading aspect ratio: 

= 1 - high aspect-ratio blading 
= 2 - medium aspect-ratio blading (default) 

= 3 - low aspect-ratio blading 

IMID indicates meanline shape: 

= 0 - meanline linear from stage 1 to stage N (default) 

= 1 - meanline constant from stage 1 to stage NMID; then linear to stage N 
= 2 - meanline linear from stage 1 to stage NMID; then constant to stage N 

The sample input file shown in table I contains four cases, each illustrating one of the new 
capabilities added to the code. Each case begins with a title card, which is indented for 
demarcation purposes. The first case is a two-stage cooled turbine with the coolant flow 
ratios, coolant temperature, and stage-efficiency specific losses included as input. The next 
three cases are a five-stage turbine. For the first of these, stage reaction (REACT) is 
specified as input (option IVD=6). The next case uses stage work factor (WF) as an input 
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(option IDIAM=4) along with exit/inlet diameter ratio to determine the stage diameters. These 
inputs were determined from the solution to the previous case and, therefore, should provide 
the same solution. The last case illustrates the use of stage work fraction (XI) as input. The 
output corresponding to this sample input is described in the following section. 

Output 

Table II presents the output that corresponds to the sample input of table I. A program 
identification title is automatically printed as the top line of the page for each new case. That 
is followed by the input title record message. The next four lines for each case are the input 
variables and their associated values. The input variable names are spelled out. The units 
for the input variable values are as described in the "Input" section. The input diameters for 
the first case are mean diameters as indicated by the MN in the variable name. Hub and tip 
diameters would be indicated by HB and TP, respectively. If diameters are calculated from 
an input work factor, as for the last two cases, the letters WF are used in the variable name. 
Note that the input diameters for the last two cases are the relative values. For a cooled 
turbine, such as the first case, two additional lines are printed to echo the coolant input 
parameters. Where stage work split is input, as for the last case, an additional line is printed 
to show the stage work fractions. 

The next group of nine lines is the computation results satisfying the input requirements. 
The output parameters are spelled out and are self-explanatory. These temperatures, 
pressures, velocities and angles are meanline values, and the continuity and efficiency 
calculations are based on these values. Note that identical solutions are achieved by taking 
the work factor and exit/inlet diameter ratio from the second case and using them as input for 
the third case to compute the absolute values for the diameters. For the cooled turbine, note 
the differences between total efficiency, which is the thermodynamic efficiency (eqn. (2)), 
rotor primary efficiency (eqn. (5)), and uncooled efficiency. These values are all the same for 
the uncooled turbine. 

The next group of four lines is the hub and tip free-vortex values of Mach number and 
angles for the last stage, where the radial variations are the largest. These flow parameters 
do not enter into the continuity and efficiency calculations, but are shown only for information. 
Following this is the meanline slope for the stages where the diameter is varying. 

Where the stage work split is input, the stage work factors are no longer assumed equal 
and the stage velocity diagrams are not necessarily geometrically similar. For this option, the 
velocity diagram parameters and efficiency for each stage are printed, as shown for the last 
case. Note the variation in stage work factor from 3.24 at the inlet diameter to 2.33 at the exit 
diameter for this case of equal stage work. 

The final output for each case are the blading geometries. Given for each stage are the 
chords, solidities, stagger angles and blade count for the stator and the rotor. Also shown for 
the last stage is the rotor blade centrifugal stress parameter AN 2 , where A is the exit annulus 
area and N is the rotative speed. 



SUMMARY OF RESULTS 


This report presents the latest modifications made to the computer code TURBAN, which 
is a preliminary sizing analysis for axial-flow turbines. The TURBAN analysis is based on 
mean-diameter flow characteristics. Program input includes flow, speed, and power or 
pressure ratio. The output presents annulus dimensions, diagram velocities and angles, and 
efficiencies. Options are provided for varying stage number, mean diameter, reaction, 
loading, diagram type, and/or work split. 

Modifications were recently made to TURBAN to satisfy user needs and convenience. 
Turbine cooling now can be accounted for in the overall energy balance and efficiency 
estimate by inputting: 

(1) ratios of cooling flow to turbine-inlet flow for each blade row; 

(2) ratio of coolant temperature to turbine-inlet temperature; and 

(3) stage-efficiency decrements due to stator and rotor cooling. 

Both thermodynamic efficiency and rotor primary efficiency are computed. 

Alternative input options have been added for defining the velocity diagrams: 

(1) specifying stage reaction to calculate stage swirl split; 

(2) specifying stage loading to calculate mean diameter; and 

(3) arbitrarily specifying stage work split. 

These options can be used in any combination. 

The Reynolds number dependency for the loss model was weakened, and an internal 
calculation of air viscosity was added as a default for convenience. The analytical modeling 
for all these modications was presented herein. 

This report also serves as an updated users manual for the modified TURBAN code. 
Program input and output are described, and sample cases illustrating the new capabilities 
are included. 
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TABLE I - SAMPLE INPUT 


COOLED HPT - 2 STAGES - Cooling parameters are input. 

& INPUT TTIN=1277 . ,PTIN=50. ,MU=2 . 69E-5 , R=53 . 37 ,GAM=1 . 38 , DIN=27 . 14, 

DEX=27 . 25 ,RPM=8285 . , POW=3630 . 7 ,W=25 . 7 , ALPHA=7 1 . 6 ,NMIN=2 ,NMAX=2 , E=1 . 17 , 
IALPH=0 , IDIAM=2 , IVD=2 , IEV=0 , IPR=0 , IU=2 , KLOSS= . 30 , ITIT=1 , 

WCOWP= .094, .070, .025,0. 0 ,TC0TP= . 48 ,DELS= . 45 ,DELR=0 . 9 SEND 
LPT - 5 STAGES - Reaction is input (IVD=6 option) 

& INPUT TTIN=750 . , PTIN=45 . , MU=1 . 55E-5 , R=53 . 37 , GAM=1 . 4 , DIN=2 1.2, DEX-25 . 0 , 
RPM=3208 . 7 , POW=5045 . 8 ,W=62 . 58 , ALPHA=6 1 . ,NMIN=5 ,NMAX=5 ,E=1 . 0 , REACT- 0 . 5 , 
IALPH=0 , IDIAM=2 , IVD=6 , IEV=0 , IPR=0 , IU=2 , IMID=2 ,NMID=3 ,KLOSS= . 30 , ITIT-1 , 


WCOWP=10*0.0 &END 

LPT - 5 STAGES - Stage work factor and relative diameters are input 
& INPUT DIN=1.0,DEX=1. 17925, IDIAM=4,WF=2. 5465, ITIT=1 &END 
LPT - 5 STAGES - Stage work split is input. 

& INPUT XI=5*0 . 2 &END 


( IDIAM=4 ) 
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TURBINE VELOCITY DIAGRAM ANALYSIS 
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